Differential Transformations of Parabolic 
Second-Order Operators in the Plane 



S.P. Tsarev * 
Institute of Mathematics, 
Siberian Federal University, 

Svobodnyi avenue, 79 

660041, Krasnoyarsk, 
e-mail: sptsarev(§mail .ru 

E. Shemyakova ^ 
Research Institute for Symbolic Computation, 
J. Kepler University, 
Altenberger Str. 69, 
Linz, Austria, 
e-mail: kathOrisc.uni-linz.ac.at 

December 17, 2008 



To Sergey Petrovich Novikov, as a development of one of his ideas. 



1 Introduction 

The theory of transformations for hyperbolic second-order equations in the plane, 
developed by Darboux, Laplace and Moutard, has many applications in classical 
differential geometry [121 IS], and beyond it in the theory of integrable systems |IH 
fT9] . These results, which were obtained for the hnear case, can be applied to non- 
linear Darboux-integrable equations [21 [TJ [151 IB]- In the last decade, numerous 
generalizations of the classical theory have been developed. Among them there 
are generalizations to the case of systems of hyperbolic equations in the plane [31 
[51 [HI 122], and generalizations to the case of hyperbolic equations with more than 
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two independent variables [9^, '23]. The non-hyperbolic case has been much less 
investigated [IHl EOl HlJ • 

Here, Darboux's classical results about transformations with differential substi- 
tutions for hyperbohc equations are extended to the case of parabolic equations. 
Thus, consider for an arbitrary solution u of the equation 

Lu = 0, L = Dl + a{x,y)D^ + b{x,y)Dy + c{x,y), b{x,y)j^O, (1) 

some Linear Partial Differential Operator (LPDO) M and a new function v{x,y) = 
Mu. One can easily compute that in the generic case v satisfies an overdetermined 
system of linear differential equations. However, there is some choice of M which 
leads to only one equation for f , namely, Liv = 0, where Li is an operator of the 
same form ([T]) allbeit with possibly different coefficients ai{x,y), Ci{x,y), hi = b. 
In this case we say that we have a differential transformation of operator L into 
operator Li with M, and denote this fact as L -—>■ Li. Also it is easy to notice that 
in this case there must exist an operator Mi such that the following equality holds: 

Ml o L = Li o M , (2) 

that is the both parts of ([2]) define the left least common multiple ILCM{L, M) in 
the ring K[D] = K[D^, Dy] of LPDOs in the plane. 
For the case of hyperbolic operators of the form 

Lh = D^Dy + a{x, y)D^ + y)Dy + c{x, y) (3) 

there are quite complete results on the possible form of the operators M that satisfy 
(EJ (see [25l Ch. VIII]): in the generic case the operator M can be determined (up to 
an arbitrary multiplier) from Mzi = 0, z = 1, . . . , /c, where Zi{x, y) are independent 
solutions of LnZi = 0. There are also some degenerate cases. As was discovered by 
Darboux, one of those degenerate cases is the classical Laplace transformation, which 
is defined by the coefficients of operator only. Relation for the "intertwining 
operator" M is widely used in the study of integrability problems in two- and one- 
dimensional cases [H m] . 

In this paper, we prove general Theorem 13.11 that provides a way to determine 
transformations L Li for parabolic equations ([T]). It turned out (Theorem 14.21) 
that transforming operators M of some higher order can be always represented as 
a composition of some first-order operators that consecutively define a series of 
transformations of the operators of the form ([T]). 

Unlike the classical case of the Laplace and Moutard transformations, the trans- 
formations considered in this paper are not invertible. In this respect the problem in 
question is analogous to the generic case that was considered in [251 Ch. VIII]) for 
operators As follows from Theorems 13. 14.21 for parabolic operators ([T]) there 
are no degenerate cases like Laplace transformations for arbitrary operators i^: any 
differential transformation of the operator ([T]) can be determined by an operator M 
of the form (ITTi) . It is of interest to consider the problem of the existence of an 
inverse transformation Li L. The order of the inverse may be higher than the 
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order of the initial transformation L — > Li. Examples show that the existence of 
such an inverse implies some differential constrains on the coefficients of the initial 
operator L. In Sec. [5] we show that these relations can imply famous integrable 
equations, in particular, the Boussinesq equation. This result is an analogue of re- 
sults [lOlITHElj for periodic chains of Laplace transformations for the operators (I3l), 
which also lead to integrable non-linear equations. 

Authors are thankful to M.V. Pavlov for useful discussions. 

2 Basic Definitions and Auxiliary Results 

Consider a field K of characteristic zero with commuting derivations dx,dy, and the 
ring of linear differential operators K[D] = K[D^, Dy], where Dj.,Dy correspond 
to the derivations dx,dy, respectively. In K[D] the variables Dx,Dy commute with 
each other, but not with elements of K. For a ^ K we have Dia = aDi + di{a). 
Any operator L G K[D\ has the form L = Ylt+j=o'^ij^x^yy "where aij G K. The 
polynomial Sym^^ = X]j+j=(i O'ijX'^Y^ in formal variables X, Y is called the (principal) 
symbol of L. 

Below we assume that the field K is differentially closed unless stated otherwise, 
that is it contains solutions of (non-linear in the generic case) differential equations 
with coefficients from K. 

Let K* denote the set of invertible elements in K. For L G K[D] and every 
g G K* consider the gauge transformation L ^ = g^^ o L o g. Then an algebraic 
differential expression I in the coefficients of L is (differential) invariant under the 
gauge transformations (we consider only these in the present paper) if it is unaltered 
by these transformations. Trivial examples of invariants are the coefficients of the 
symbol of an operator. A generating set of invariants is a set using which all possible 
differential invariants can be expressed. 

Theorem 2.1. [T7[ \2BI The action of the gauge group on operators of the form ([T]) 
has the following generating system of invariants: 

h = b, 

h = Cx - aax/2 - bay/2 - axx/2 
+ {bxa^/A-bxC + bxax/2)/b . 

Note that if an operator ([T]) has only constant coefficients then Ji is a constant 
and I2 = 0. If the field of coefficients K contains quadratures (differentially closed), 
it is easy to prove the inverse statement: 

Proposition 2.2. Let the field of coefficients K be differentially closed. The equiv- 
alence class of ([T]) with respect to gauge transformations contains an operator with 
constant coefficients if and only if Ii is a constant and I2 = 0. 

Proof. Let Ii = b have a constant value and I2 = 0. Consider an operator L = 
Dl + aDx + bDy + c from the equivalence class. Using the gauge transformation with 
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g = exp ( — I / o dx^ one can make a = 0. Then I2 = implies = Cj; — hxc/h. Since 
Ji = 6 is a constant, we have c = c{y). Applying the gauge transformation with 
g = e-^"^/^*^^ to L we obtain = + bDy, which has constant coefficients. □ 
So every operator ([T]) with constant Ii = b and I2 = can be transformed into 
operator D"^ + Dy using substitution ?/ i-^ const ■ y and gauge transformations. 

Lemma 2.3. Without loss of generality one can divide the symbols Sym(M) = 
Sym(Mi) by any non-zero g E K . The operator L and the symbol of Li are left 
unchanged. 

Proof. Indeed, multiply the both sides of (E]) by 1/g on the left: o L = 

^LiQ o iMi = Lf o iMi. Then "new" M and Mi have the coefficients of the "old" 
ones divided by g, while Li is subjected to the gauge transformation with g, and, 
therefore, its symbol is unchanged, while the other coefficients can be changed. □ 

Lemma 2.4 (Simplification by gauge transformations). In one can assume 
without loss of generality that a = 0, that is there exists a gauge transformation 
that transforms L, M , Li and Mi into operators of the same form such that the 
coefficient of L at is 0, and the equality M o L = Li o Mi (|2]) is preserved. 

Proof. It is enough to apply the gauge transformation with g = exp( — | f adx) 
to all operators in ([2]). This gauge transformation do not alter the symbols of the 
operators, and, therefore, does not interfere with the simplifications from Lemma [^7S1 

□ 



3 First-Order Transformations 

Consider L of the form ([T]) and an operator Li of the same form: Li = D"^ + 
ai{x,y)Dx + bi{x,y)Dy + ci{x,y). Then a differential transformation of the ffist- 
order that transforms L into Li exists if there exist 

M = p{x,y)D^ + q{x,y)Dy + r{x,y) , 
Ml = pi{x,y)D.^ + qi{x,y)Dy + ri{x,y) 

such that ([21) holds. The comparison of the symbols implies pi = p, qi = q- 
First consider the case p 7^ 0, q 7^ . 

By lemma 12.31 without loss of generality one can assume p = 1, and a = 
by lemma 12. 4[ Equating the coefficients in ([2]) we have ai = — 2y , 61 = 6, ci = 

{—2bqx + bxq + q'^c-\-q%y-\-2q^ — bqyq — qxx(l)/(f'^ fi = r — 2(lng)3;, and two constrains 
on the coefficients of the operators L and M: Ci = 0, Cq = 0, where 

Co = -2gcg^ + c^q^ + q^Cy + 2r6g^. - rb^q - rq\ - 2rql + 

+rbqyq + rq^^q + 2qxrxq - bvyq^ - r^^q'^ , (4) 
Ci = -2bqx + b^q + q%y + 2ql - bqyq - q^^q - 2g^rg + 2r^g^ . (5) 
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We see from (j4]), ([5]) that given the coefficients of the operator L, one can always 
find solutions r, g of these equations in the differentially closed field K, that is ev- 
ery operator ([T]) admits infinitely many transformations with different operators M. 
The equations (jlj), ([5]) for r, q can be solved explicitly with the help of two arbi- 
trary (independent) generic solutions of the equation ([T]). Indeed, given a ffist-order 
operator M that satisfies the constrain ([2]), the following system of equations 



is consistent and has a two-dimensional space of solutions, which is parameterized, 
for example, by the values u{xo, yo), Uy{xo, yo). In fact, we can express the derivatives 
of u of any order with respect to x in terms of its derivatives with respect to y from 
the second equation Mu = 0. Substituting those into the ffist equation Lu = 
0, we have an expression for the second derivative Uyy, provided q ^ 0. On the 
other hand the consistency of is guaranteed by ([2]), which can be rewritten as 
qDyLu — D^Mu = mod (L, M). Conversely, a basis zi{x, y), Z2{x, y) in the space 
of solutions of (EI) allows us to reconstruct M: the conditions Mzi = 0, Mz2 = 
give a system of two linear algebraic equations for the coefficients r, q, and we can 
easily determine the operator M: 



Since the values Zi{xo,yo), {zi)y{xo,yo) are lineally independent, the denominator of 
this expression is non-zero. 

Vice versa, the choice of two arbitrary lineally independent solutions Zi, z^ of 
the equation ([T]) defines the operator M by the formula ([7]). The operator M in 
its turn implies a differential transformation of L, that is the equality Indeed, 
compute the derivatives fy, Vxx of the function v = Mu for an arbitrary solution 
u of the equation ([T]), then using ([T]) we can remove all the terms that contain Uxx, 
Uxxx, Uxxy Using an appropriate combination Lv = Vxx + o.i{x,y)vx + bi{x,y)vy 
we can also remove the terms with Uxy, Uyy, leaving Ux, Uy, u only. Since the 
expression Lv vanishes after the substitution u = Zi it must be proportional to Mu: 
Lv = L{Mu) = ci{x,y)Mu, which implies ([2]) with Li = L — ci for an arbitrary 
function u{x, y). 

Note that in the considered case the coefficients at Dy, Dx in M are non-zero. 
From now on we refer to such transformations as X -|- gF-transformations. Below 
we consider the cases when one or another of the coefficients is zero separately. 
Therefore, we will prove the following statement: 

Theorem 3.1. For every operator L = + aDy + bDy + c there exist infinitely 
many differential transformations with operators M = Dx + q{x,y)Dy + r{x,y). If 
q then the operator M is defined by the conditions Mzi = 0, Mz2 = 0, where Zi 
are two arbitrary chosen independent solutions of the equation (Qp. The operators 




(6) 



Mu 



{Zl)y {Zl)x 

^2 (^2)5/ (^2)0: 



Zl {Zl)y 

Z2 {Z2)y 



(7) 
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of the form M = + r{x,y) are defined by the choice of one solution Z\ of the 
equation (QP and by the condition Mzi = 0. The intertwining operator of the form 
M = Dy + r(x, y) does not exist for generic L. 

The degenerate cases of operators M of forms M = + r and M = Dy + r are 
considered below. 

Case p^O, q = (M = D^ + r) 

Without loss of generality one can assume p = 1 and a = 0. If we equate the 
corresponding coefficients in ([2]), we have ai = — lia{b)x, bi = b, ci = c+r\n{b)x — '2rx, 
ri = r — ln(6)^. and an equation 

= -c ln(6)^. + c^. - ln(6)^ + 2rr^ + ln(6)^. - bvy - r^^ , (8) 

for r. We apply the same trick as in the non-degenerate case in order to determine 
the operator M in terms of solutions of the initial equation ([T]). Now we choose one 
solution Zl and require M to satisfy the condition Mzi = 0. We get 



M{u) 



U Ux 



(9) 



Indeed, given an operator M such that the intertwining equality (Ej) holds, an ap- 
propriate Zl is found as a solution of the consistent system ([6]), which now has a 
one-dimensional solution space. 

Conversely, given a solution zi of the equation ([T]), M can be found from Q, 
then for v = Mu the derivatives Vx, Vy, v^x are simplified using ([T]). 

Then an appropriate combination Lv = Vxx + (ii{x, y)vx + bi{x, y)vy contains only 
Ux and u (there are no terms with Uyy\). The obtained expression Lv vanishes if we 
substitute u = Zi and therefore it must be proportional to Mu, which implies 

Later on we refer to such transformations as ^-transformations. 

Case ^ = 0, ci^Q (M = Dy + r) 

Without loss of generality we can assume g = 1, a = 0. If we equate the 
corresponding coefficients in ([21), we obtain in particular Tx = 0, Cy — rby — bvy = 0. 
Thus, r = r{y) can be found only for some particular functions b, c and for an 
arbitrarily chosen L = + aDx + bDy + c there is no differential transformations 
with M = Dy + r. 

Notice also that an attempt to construct M by the formula 



M{u) 



U Uy 

Zl izi), 



z7\ 



would not lead to any success either: for such an operator M and v = Mu the 
derivatives Vx, Vy, Vxx simplified with ([T]) would contain Uxy, Uyy, Ux, Uy, u, and we 
cannot not find an appropriate combination Lv = Vxx + cii{x, y)vx + bi{x, y)vy having 
only Uy, u. 

Therefore, Theorem 13.11 is proved. 
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Note that when differential transformations with M = + qDy + r are apphed 
to the operator ([T]), the new values of the basic invariants (that is the values of 
invariants Ji and I2 for Li) are 

II = Ii = b, 

II = h- 2bq,Jq^ - bl/{qb) - b,by/b + b,Jq + b,y - b.qjq^ + Aqlb/q^ . 

When differential transformations with M = + r are applied the new values 
of the basic invariants are 

II = Ii = b, 

II = l2-l/A{8b''r,,-12b,r,b^ + 8rbbl 

-4rbX^ + 2b%yb, - 9bl - 2b%,,^ - m^bb,, - 2b%y)/b^ . 

Example 3.2. Consider an operator 

L-D +^^±^D -A 

The equation L{z) = has the following solutions Zi = x'^, Z2 = x + y. Using the 
determinantal formula ([7]) compute 

X X 

and M^ = D, + 



x{x + 2y) X x"^ ^ {x + 2y)x {x + 2y)x'^ 

Note that Li cannot be obtained from L by any gauge transformation. Indeed, 
the value of the invariant I2 for L is /2 = -27^— r, while the value of I2 for Li is 

X \x-\-y) 

jl ^ 2(a:^-2xy-V) 
2 x{x+y){x+2y)^ 

Example 3.3. Applying the differential transformation with M = D^ + q{x, y)Dy + 
r{x, y) to L = D^ + Dy (provided conditions (HI) and ([5]) are satisfied or equivalently, 
provided M is in the form (jTj)) we have Mi = + q{x, y)Dy + r — 2(lng)a; and 

L = Dl + Dy > Li = Dl- 2qJqD.^ + Dy + (qyq + q-^-^q - 2ql + 2q^)/q^ 

l2 = ll = -2q,,/q^ + Aql/q' . 

Example 3.4. Applying the differential transformation with M = + r{x,y) to 
L = Dl + Dy (provided the condition (IHl) is satisfied or equivalently, provided M is 
in the form([n])) we have Mi = M and 

L = Dl + Dy ^ Li = Dl + Dy-2r, , 

/2 = ^ 4 = -2r,, . 
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4 Transformations of Arbitrary Order 



We show that differential transformations of arbitrary order of a generic operator ([T]) 
can be expressed in terms of some number of partial solutions of ([T]). In [25; Ch.VIII] 
analogous formulae were introduced for hyperbolic operators (I3l). 

First of all, given some transforming operator M of higher order satisfying 
we can use the operator L to remove all terms having derivatives with respect to 
y (generally speaking, this manipulation increases the order of M). The resulting 
operator has the form 



M 



i=0 



(10) 



Below we call the corresponding transformation an {m) -transformation. 



Theorem 4.1. Given an operator ([T]) and m lineally independent generic partial 
solutions zi,...,Zm of the corresponding equation L{z) = 0, then there exists a 
differential transformation with 



Mu = {}{x, y) 



du 






dx"^ 




d"'zi 


dx 


dx"^ 


dZm 


d^Zm 


dx 


dx^ 



where d{x,y) is arbitrary. Conversely, every {m) -transformation of an operator of 
the form ([T]) corresponds to some operator M of the form fill I) . 



Proof. Having computed the derivatives Vx-, Vy, Vxx oi v = Mu for an arbitrary 
solution u of equation ([1]), we use ([1]) as above to remove all terms that contain 
derivatives with respect to y. The remaining terms will contain only some linear 
combinations of the derivatives -D^m, s = 0, . . . , m + 2. Choosing some appropriate 
combination Lv = Vxx + 0'i{x,y)vx + bi{x,y)vy we can remove terms with 
D!^^^u, and leave terms with D^u, s = 0, . . . ,m only. Since the resulting expression 
Lv vanishes when we substitute any u = Zi, we conclude that it must be proportional 
to Mu: Lv = L{Mu) = Ci{x,y)Mu, which implies ([2]) with Li = L — ci for 
an arbitrary function u{x,y). The only requirement is the non- vanishing of the 
Wronskian det{Dizi), i = 1, . . . ,m, j = 0, . . . ,m — 1. 

Conversely, given the intertwining operator M of the form ( fTOl) satisfying ([21), 
consider the system ([6]) . The consistency of the system is equivalent to ([2]) , which al- 
lows us to choose a basis of its m solutions with non- vanishing Wronskian det{Dlzi), 
i = 1, . . . ,m, j = 0, . . . , m — 1, and obtain the required form (fTTj) of the operator 
M. □ 

Theorem 4.2. An arbitrary (m) -transformation of an operator (J\) with m > 1 can 
be represented as a composition of first- order differential transformations. 



8 



Proof. Consider an operator M in the form (ITT!) and the corresponding solutions Zj. 
Then Zi generates a first-order transformation with M of the form © , which trans- 
forms L into some L of the same form ([T]). Others Zi, i = 2, . . . ,m are transformed 
into solutions Zi = Mzi of the equation Lz = 0. Since Mzi = 0, Mzi = 0, then if we 
divide the ordinary differential operator M by M, the remainder is zero: M = PM, 
P e K[Dx]. implies that the operator ILCM{L,M) = MiL = LiM is divisible 
by ILCM{L, M) = M^L = LM, that is M^L = N^MiL = NiLM = L^M = L^PM, 

which implies NiL = LiP. Thus we have obtained an intertwining operator P, 

p 

whose order is less by one, such that L — > Li. The induction by the order m of 
the intertwining operator completes the proof. □ 

5 Generalized Moutard Transformations and Dif- 
ferential Transformations. Periodical Differen- 
tial Transformations 

An important subclass of the considered class of the parabolic operators are opera- 
tors 

L = Dl-Dy + c{x,y) . (12) 

In [H], a modification of Moutard transformations for such operators was 
suggested and applications to the construction of solutions in the Kadomtsev — 
Petviashvili (KP) hierarchy of equations were given. As we show below, some of the 
examples considered in [8J can also be obtained by our method. Direct application 
of the above results proves the following lemma. 

Lemma 5.1. X -transformations preserve the class of the operators ( fT2l) . For M = 
Dx + r{x,y) the condition (IHl) for the existence of such transformations has the 
following form: 

Cx + 2rrx + ry-rxx = , (13) 

and 

Mi = M , Li = Dl- Dy + c-2rx . (14) 

The basic invariant I2 transforms as follows: I2 = Cj. — > J| = c^. — 2r^a.. // the 
operator M is given in the form ^ for some partial solution Zi = Zi{x,y) of the 
equation L{z) = 0, we have 

r _ n2 n '^^ix " ^l^ly " ^\xxZ\ 

^1 - ^x - j2 • 

^1 

Note that X + gF-transformations do not preserve the class of operators ( JT2ll : 

Example 5.2. The equation {Dl — Dy)z = has partial solutions Zi = x, Z2 = e^^^ . 
The formula ^ implies M = + (1 - l) - i and L, = Dl- ^^D, - Dy - 
^(^^^^-^ ■ However, the gauge transformation with g — {x — l)/x reduces Li to the 
form dni): Lf = Dl - D 



y 
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This example and the one below show that classical examples of functions c{x, y) 
obtained in [B] can also be obtained by the application of one or several differential 
transformations. Actually, both approaches can be considered as two-dimensional 
generalizations of Darboux transformations for the one-dimensional Schrodinger op- 
erator — c{x, y). 

Example 5.3. Consider a differential transformation of L = + Dy with M = 
Dx + r{x,y). Choosing r = 1/2 — tanh(a; -|- y) satisfying the condition (fT3l) of the 
existence of the transformation, we have 

L = Dl + Dy L, = Dl + Dy+ ^ 



cosh(a; -|- y)"^ 
. 4sinh(x + ?/) 

i2 — U > In — TT TT • 

cosh(a; + yy 

Now we study the invertibility of a given transformation L Li, that is, the 
possibility of finding a transformation Li L, possibly of higher order. 

Example 5.4. X-transformation of the operator L = — Dy — + 2x with 
M = Dx + x"^ results in the following operator: Li = — Dy — x^ — 2x. This 
transformation has the inverse X-transformation with N = Dx — x^. 

As the simplest examples show, an inverse transformation does not exist for a 
generic operator L. In fact the existence of an inverse transformation implies a 
system of constrains on the coefficients of L. In some cases, it produces known 
integrable equations. First, Theorem 14.21 implies that the existence of an inverse 
transformation, that is the existence of a composition L ^^-^ L, is equivalent to the 
existence of a transformation P = N-M of higher order that transforms the operator 
L into itself: Pi ■ L = L ■ P. For operators ffT^ the existence of such an operator 
implies a particular case of the standard problem of classification of Lax pairs: for 
P of order one or two of the form ( |T0|) this leads to potentials c(x, y) of simple form; 
the existence of an operator P = ps^x, y)D^ + P2{x, y)Dl + y)Dx + Po{x, y) of 
the third order implies Pi = P and P = AD^ + 6cDx + Po{x, y) (up to some simple 
transformations) and the system 

{Pq)x = "iiCy -\- Cxx) , ^-^^-^ 
iPo) y '^Cxy GcCx Cxxxj 

that is the famous Boussinesq equation for c: 

Cyy (C -|- Cxx /^)xx- 

The system ( 1T5|) coincides with the well-known representation ([H formula (7)]) for 
the Kadomtsev-Petviashvili equation in the stationary case t/t = 0, which gives the 
Boussinesq equation. 
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